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ABSTRACT 


In this thesis a model of a 12 pulse converter is 
described for dynamic digital simulation of HVDC system. 

Both constant current and constant extinction angle controls 
are represented in detail, 

A computer program has been developed based on these 
models, for the digital simulation of a two terminal HVDC 
system with six or twelve pulse converters at both ends. 

In addition to the converter and its controls DC and AC 
system along with filters are also represented in detail. 

Test simulations are carried out both for steady state 
and transient conditions to illustrate controller 
performance. 



CHAPTER 1 


INTRODUCTION 

1.1 GENERAL 

Despite the general acceptance of ac transmission, dc 
transmission made its successful beginning in 1954 with 
the Gotland link (Sweden), leading to an awakening of 
interest in HVDC transmission. The advent of thyristors 
and other solid state devices made converter control simple 
and reliable thus, significantly adding to the advantages 
of dc transmission. Hence, in recent years, HVDC transmission | 
has become a viable alternative over ac transmission for | 

i) bulk pow::r long distance transmission and ii) inter- 
connection of ac systems. 

To augment the existing power system with dc link, 
the study of operation and control of HVDC systems for planning j 

■ I 

and design of composite ac/dc system becomes necessary. These I 

i 

studies are usually conducted on HVDC simulator. Hov^ever, i 

the HVDC simulator has less flexibility and involves increased I 

i 

cost in representing wider range of system configuration and 
control schemes. This has given rise to the development of 
analytical techniques and digital simulation. The extent of 
mathematical modelling of HVDC system depends on the nature 
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of study. Before describing the objectives and outline of 
this thesis a brief review of the literature pertaining to 
digital simulation of HVDC systems is presented, 

1.2 REVIEW 

Dynamic digital simulation involves the study of system 
responses over a vdde range of frequency ranging from 
fundamental to higher, both under normal and abnormal condi- 
tions. Since the converter has a time varying topology 
(i,e, the conduction pattern changing with, the firing of the 
valve and its cessation), the major difficulty lies in 
representing the converter to suit the nature of studies. 

There are two approaches to the problem of handling 
the time varying topology of the convertor. 

The first method is based on the use of tensor analysis 
technique [l]. Based on the conduction pattern, a transfor- 
mation matrix is generated by the digital computer, using 
which converter dynamic equations corresponding to each 
state of the converter are formulated. Hingorani et ai. [2] 
represented a two terminal dc link based on this approach. 
Later, Hay et al, [3] introduced a general method for 
representing a multi converter system for two terminal dc 
link, Williams et al, [4] have formed the loop impedance 
matrix of the converter network assuming all the valves 
of the bridge to be conducting. Depending on the conduction 
pattern of the converter the loop impedance matrix is 
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reduced to get the connection matrix. Vovus et al, [5] 
extended the above method to simulate 12 pulse 2 terminal 
HVDC system, Padiyar et al, [6] represented the converter 
based on graph theoretic approach, assuming the dc link 
current to be continuous. The topological analysis leads 
directly to the derivation of converter dynamic equations 
and express ions for dependent variables such as valve 
voltages. 

The second method is based on the storing the sets of 
differential equations describing the converter topology 
for different valve conduction patterns. Corresponding to 
the conduction state, the appropriate set of equations is 
chosen at any given time. El~Serafi et al, [7] formulated 
a set of five general equations in phase quantities for the 
normal modes of operation of a converter. The equation for 
the particular mode of operation is obtained by substituting 
certain empirical conditions in the general equations. In 
the central process method of Hingorani et al, [8], a central 
process starts from the instant of firing a valve and 
continues upto the instant of next firing. The continuous 
converter operation is simulated by successive use of the 
central process which is a set of differential and boolean 
equations. Depending upon the state of converter the 
appropriate differential equations are selected and solved , 
The method has been used to simulate a two terminal dc link 
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assuming normal converter operation. Htsui et al, [9] 
modelled a converter by formulating a set of generalised 
differential equations which are derived corresponding to 
all the possible normal combinations of conducting valves. 

There is an alternative approach whereby the converter 
is represented by a network of fixed topology but time 
varying impedances [lO], 

A valve is represented by high impedance when it is off 
and low impedance when it is on. Though this method is‘ 
conceptually simple, it needs a small integration step size 
for the solution of the converter state equations due to 
time dependent impedances. 

The design of control system is one of the aims of 
dynamic digital simulation. Very few papers have dealt with 
this aspect in detail. References [2,5,8] algebraically 
relate firing angle to the current error. This simplification 
may not be adequate when ac system is represented in detail. 
Reeve et al, [H] modelled constant current controller and 
extinction angle controller based on pulse phase control 
firing scheme and illustrated their capability by simulating 
a multiterminal HVDC system under transient conditions, 

Padiyar et al, [6] represented the constant current controller 
using individual phase control (IPC) or EPC (using pulse 
frequency control) firing schemes and the constant extinction 
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angle controller using IPC scheme. They validated the 
models by performing transient studies on a raultiterminal 
HVDC system. 

Vovas et al, [5] considered the detailed representa- 
tion of ac system and simulated ac line faults on a parallel 
'«c/dc. link connecting two areas. Each area includes 
detail representation of synchronous machine, El-Serafi [7] 
investigated the transient behaviour of a synchronous 
machine connected to a converter. Hay et al, [3j represented 
the ac system as a simple impedance calculated from the 
short circuit I'atio (SCR) at the converter bus and the 
impedance angle of the ac system. In the references [2,5,7] 
the harmonic filters have also been included in simulation, 
the study of harmonic interaction between the ac and dc 
systems is reported in the papers [12,13], 

1.3 THE OBJECTIVE OF THE THESIS 

The objectives of the thesis axe 

i) to report a methodology of representing a HVDC converter 
based on equivalent circuit approach using graph theore- 
tic analysis. This representation can be used for 6 as 
well as 12 pulse converter | 
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ii) to develop a digital computer program to simulate 

a two terminal dc link based on the proposed converter 
model and v;ith detailed representation of such 
as converter controls, DC and AC system alongwith filters? 

iii) to simulate a two terminal HVDC system under both 
steady state and transient conditions and to observe 
the effect of IPC and EPC firing schemes on the system 
response using the computer program, 

1.4 OUTLINE OF THE THESIS 

The chapterwise summary of the thesis is as follows i 

i) Chapter 2 deals with model formulation of dc sub- 
system consisting of the converter, dc line, current 
controller and the extinction angle controller. 

ii) Description of the development of the computer program 
alongwith its flow-chart is given in the Chapter 3, 
Steady state and transient performance of a HYDC link 
without detailed ac system representation have been 
presented, 

iii) The state space model of the ac system and harmonic 
filters is developed in Chapter 4, The converter is 
treated as a dependent current source for the ac system 
whose effect is represented as a dependent voltage 
source in the convertor model. In interfacing the ac 
and dc side variables the identity of transformer 
connection and tap is maintained. The computer program 
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developed in the Chapter 3 is augmented to include 
the ac system model and its capability is validated 
using single and two terminal HVDC system. The effect 
of harmonic filters and effectiveness of IPC and EPC 
firing schemes on the system performance are studied, 
iv) Chapter • 5 , reviews the major contribution of the 
thasis and gives the future scope of vrork. 
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CHAPTER 2 

MATHEMATICAL MODELLING OF A 12-PULSE CONVERTER 

2.1 INTRODUCTION 

A HVDC system mainly consists of (i) the ac system 
feeding the converter (ii) converters with their associated 
controls and (iii) the DC network consisting of transmission 
lines. Each component is modelled separately and in a modular 
fashion. The individual component models are interfaced 
using appropriate variables. 

This chapter deals with the 12-pulse converter repre- 
sentation alongv/ith its associated control system. The 
converter is modelled as a variable voltage source behind 
a variable impedance. The formulation of the converter 
equations corresponding to all the possible modes, with the 
assumption of continuous dc link current, becomes simple 
using graph theoretic approach. The converter equations are 
formulated using cut set matrix and the constituent equations 
for the elements of the associated graph. The dynamics of 
the converter control using both individual phase control 
and equidistant pulse control firing schemes, are represen- 
ted and simulated in detail, 

2.2 COWERTHR REPRESENTATION 

The modelling of the 12-pulse converter has been 
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achieved through the modelling of the individual 6-pulse 
converter unit. Since, the extension of the 6-pulse 
converter model is dealt later in this article, the formu- 
lation of 6-pulse converter model is taken up first. 

A three phase bridge converter system is shown in 
Fig, 2,1, This includes the resistances and inductances 
of the converter transformer L^,L. ,L^) and the 

dc smoothing reactor Both the ac and the dc voltage 

sources are not constant and are actually the output of ac 
and dc network models. The effect of the converter on the 
ac and dc networks are represented by the injection of 
currents into the respective network. 

The graph associated with the converter of Fig, 2,1 is 
shown in Fig, 2,2, In total there are 9 elements of which, 
the first 6 represents the values and the orientation of 
each of these is along that of the direction of current 
flow in the valve. Elements 7 and 9 correspond to the 
equivalent circuit representation [l] of the ac system 
feeding the converter. The element 8 includes the series 
combination of and V^, In deriving the graph of 

Fig, 2,2 the grading and damiping circuits across the valves 
are ignored. The tree of the graph in Fig, 2,2 will have 
4 elements. Although there is no unique tree for a graph, 
in v;hat follov>/s elements 7 and 9 are chosen as two of the 
elements of tho tree. In addition, any two of the conducting 
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valve elements are included in the tree, (It is to be 
noted that there are, always, at least two valves conducting i 
dc current is to be continuous. 

The constituent equations of elements 7 and 9 can be 
written as 


“T1 " ~11 


( 2 . 1 ) 


and that of element 8 is 




( 2 . 2 ) 


where 


<^Ti = [Vy Vg] ^ I Vj is the voltage across 

the element j 


i^^ = [i<^ is ihe current through 

element j 


^8 “■ ^d ^d^ 


Go = 

8 c 


and 


C®c-®b ®a~®bJ 


the superscript X indicates transpose and s^^^b ®c 
line to neutral ac voltages of phases a,b and c respectively, 
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Zi, = R,,(l+cp) 

Here c,p and are defined as 

L, L, 

^ = if = rT = r ’ P - 

a D c 


«ii= 


R,+Ph, 

a D 


R. 


R. 


c b 


c, the time constant of source impedance is assumed to 
be the same in all the phases* 

The elements 1 to 6 represent the valves and are 
governed by the following equations; 


V,, = 0 if the valve is conducting 

(2.3) 

ij^ = 0 if the valve is non-conducting 

Vj^ and are the voltage across and the current through the 
V al ve . 

From graph theory, the branch and link variables are related 
by the follov<;ing expressions 


^T1 


1 

00 


= -B 

\2 

: ^T2 

i 

. \3 

^ J 




(2.4) 
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hi 


r 


T 

^T1 ■ 

\2, 

= 

« - 

■ 


^T2 V 

;^L3j 

J 

r 

J 


where subscript T and L refer to branches and links 
respectively. B is the fundamental cutset matrix 
corresponding to link elements. 

The tree branches are partitioned into two sets, one 
set (T^) consisting of elements 7 and 9 and the other (T 2 ) 
consisting of the conducting valves. The link elements are 
partitioned into 3 sets. The first set (L^) corresponds 
to element 8 and the other two sets (L 2 and L^) consist of 
conducting and non-conducting valves in the links respectively. 
The matrix B can then be pax-ti tinned as 


B - 

h 

\ll J 

il2 

i ®L13 

- 4- 



■^L2l ; 

®L22 

1 ®L23 


Using (2,3) in conjunction with (2,4) and (2,5) it follows 
that 

is = ° 

v^2 = ° (2.6) 
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Expanding equations (2.4) and (2.5) and using (2.6), the 
following relations are obtained 


^T1 

“\ll^8 " 

■ \l2^L2 

(2.7) 

^T2 ~ 

"■®L21^8 ■ 

' ®L22^L2 

(2.8) 

Vr = 

T 

3^ V 

"^LirTi 


(2.9) 

\2 "" 

° = <ii2' 

r 

T1 

(2.10) 

II 

CO 

> 

®L 13^11 


(2.11) 

tion (2. 

9) can be 

modified using (2,7) 

and (2,1) as 

^8 ^ 

= -Z i 

18 

^2--L2 ^11®11 

(2.12) 


and can be written by using (2,10), (2*7) and (2,1) as 
(l+cp)ij^2 = ”^4 ^3^8 ^4 ®L12®11 (2.13) 

where 

^1 “ ®L11^iAi 1^ R^(l+cp) 

^2 = ^iAAi 2(1-^"P^ = R2(l+cp) 

^3 “ ®L12^iAi 1^ R3(l+cp) 

^4 ~ ®Ll2^1l\l2^ ~ R4(l+cp) 

Equation (2.12) can be further simplified by using (2,13) 
to get 
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V 


8 




(2.14) 


where 


^eq = (Ri-R 2 “;^R 3 )(l+cp) = r^(l+cp) 


®eq “ ®11 


Vg is the voltage across the clement 8, which is the dc 
voltage before the smoothing reactor of the bridge. Thus, 
using (2,14) an equivalent circuit (Fig, 2.3) can be formed 
where and are the impedance and voltage source 

respectively in the equivalent circuit. 


The complete model of a 6-pulse converter is thus as 
shown in Fig, 2,4, The dynamic equation which gives the dc 
current of the converter is 


fj, 

dt " 


r ^8 


(2,15) 


where 

^ “ ^eq + J L = c Tgg + 

In addition to (2.15), equation (2,13) has to be solved for 
the case of three and four valve conduction mode. The 
currents in the conducting valves in the branches and links 
are given by the equations (2,8) and (2,13) respectively, 
vdiereas the voltage across the non-conducting valves is given 
by the equation (2.11). 
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FIG. 2.4: EQUIVALENT CliiCUIT OF 6 PULSE CONVERTER 



FIG. 2.5: TFAlISFDRiVER CONNECTION FOR 12 PULSE 
OPERATION 
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Example j 

An example of 3 valve conduction mode is taken up 

to illustrate the above- analysis* Assuming that the valves 

1,2 and 3 are conducting and out of which the valves 2 and 3 

are in the branches (set T^) and the valve 1 in the link* 

The remaining valves form the non-conducting set, L^. The 

matrix B is partitioned accordingly and is shown below, 

8 1 6 5 4 

! 

I 

1 

-i - 

} 

I 
I 

1 

^1*^2’ ^3 -^4 computed from the expressions 

/ defined above. Thus, the equivalent circuit parameters are 

^eq ^1~V4^^3 * 

From the knowledge of tlie sub-matrices of the matrix B and 
knowing r„„ and , the dynamic equation giving the 

0q ' 

current through the conducting valves in the link is 
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and the dynamic equation for the dc current of the converter 

pig = - a ig + (e^<,-V^)/L 

R and L are defined earlier. 

The expressions of the currents through the conducting 
valves in the tree and the voltages across the non- 
conducting valves are 

H ** h • 

V4 m Hg^Cigtcpig) + *«r®a 

These are computed after solving the dynamic equations, 

2,2.1 Extension of the model to 12-pulse converter 

Two six pulse bridge in series constitutes a l2-pulse 
converter. Povyer to each' of the bridges are supplied from 
two transformers, one having v^ye/wye connection and the other 
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delta/ wye as shown in Fig, 2,5, The above mode of transformer 
connection is adopted in order to produce tv;o sets of voltages 
having a relative phase difference of 30 degrees between the 
coi'responding phases. Each 6 pulse bridge can be represented 
by the equivalent circuit of Fig, 2,3 and the series connection 
of two such circuits with and gives the model for 

12 pulse conve'rter shown in Fig, 2,6, The circuit in Fig, 2, 6 
can be reduced to that of Fig. 2*7, The equivalent circuit 
parameters of this reduced 12 pulse converter model are given 
by 


E 


eq 


R. 


2 


'eqj 


eq 


3=1 


■eqj 


3=1 


^3 ^eqj 


where e and are equivalent circuit parameters for the 

eq3 QQ3 


bridge j 


The primary dynamic equation of the 12 pulse converter is 


di 


8 


u’ 


dt 


P- is + 


eq 


( 2 . 16 ) 


where L’ = | R« = 




Fig, 2,6 INTERCOMNECTION OF TWO EQUIVALENT CIRCUIT 
MODELS OF 6 PULSE CONVERTER TO FORM 12 

PULSE CONVERTER 



Fig. 2.7 REDUCED EQUIVALENT CIRCUIT MODEL OF 
12 PULSE CONVERTER 
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The equation (2.16) is of the same form as (2.15) except 
for the change in the equivalent circuit parameters. At 
every instant, (2,16) has to be solved for the 12~pulse 
converter and for each 6 pulse bridge, (2,13) has to be 
solved. The matrix B associated with each bridge changes 
whenever a valve ceases or commences conduction. The new 
matrix B is computed each time depending on the tree 
elements. 

The procedure outlined above can be compared with the 
method described in [2,4,6] which are similar in objective. 
The method given here is simpler both conceptually and 
computationally. The main difference of this 6 pulse model 
with [6] is that the present model can be easily extended 
to represent a 12 pulse converter. Inclusion of a different 
valve characteristic such as finite impedance during the 
non-conducting period of the valves other than that of an 
ideal switch can also be accomplished. In doing so, the 
constituent equation (2,i) would be different. In contrast, 
the reference [2,4] implicitly assume that the valve 
impedance is infinite when it is non-conducting. The number 
of state equations per converter terminal will vary depend- 
ing on the number of conducting valves per bridge (N^) and 

b 

is given by 1 + ^ (Nj_"-2), where b is the number of 

i=l 

bridges per converter terminal. 
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2.3 CONVERTER CONTROL 

Each converter in a IT/DC system is equipped with a 
controller that determines the instant of firing of each 
valve. The firing pulses are generated from a control 
signal which is the output of the current or extinction 
angle controllers. Normally current control is used at the 
rectifier and constant extinction angle (CEA) at the inverter. 
Basically, there are two firing control schemes i 

i) Individual phase control (IPC), vdiere the firing pulse 
for each valve is determined with respect to its own 
commutation voltage. 

ii) Equidistant pulse control (EPC), Vv’here the controller 
generates pulses at equal intervals independent of 
the distortion in the ac voltage waveforms. 

2,3.1 Firing Control Scheme with IPC 

The details of this firing pulse generation scheme are 
described in reference [15] for analogue implementation. 

The same concept, for pulse generation is used here on the 
digital computer. A converter is generally equipped with 
both constant extinction angle controller and constant 
current controller in order to permit change in the mode of 
operation. 
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The control signal of the converter as derived in 
reference [i5] is 


^0 = 


wt 

/ 


e^.(x)dT + 2X^1^ + 


(2.17) 


where e . is the commutation voltage for valve 
commutating reactance and is the dc current 
output of the first order current controller, 


3> is 

• is 

cc 

expressed 


the 

the 

as 


V 


cc 


K 

l+T^p 


^^ref ” 


(2.18) 


where K is the controller gain, T^ is the controller time- 
constant and lpg£ is the current reference. The state equat 
equation of the controller can bo written from the equation 
(2,18) and solved for V , 


Owing to harmonics, the ac voltage at the converter 
bus may not be sinusoidal. Thus, for the purpose of simula- 
tion the numerical integration is employed to evaluate the 
integral in equation (2,17). However, in case of sinusoidal 
ac voltages a closed form, expression can be derived for the 
integral [l5]. 


The firing pulse is initiated at the positive going 
zero crossing of the control signal E To ensure minimum 
extinction angle operation of the converter, the voltage, V 

C w 

is clamped to zero while going negative. With increase in 

V , the angle of advance (p) increases for the converter and 
0 c 
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if is large enough, p may increase to 180°, Hence, the 
control of the converter is over the full range of operation 
(rectification and. inversion). In some extreme cases, there 
lies a possibility of the signal V becoming larger thus 
forcing the control signal E to be positive always. This is 
taken care of by adding a negative pulse of duration equal 
to minimum delay angle ( to the negative peak of the 
integral in equation (2,17), The amplitude of the pulse is 
large enough to ascertain that the control signal E goes 
negative for atleast the duration of the pulse and ensuring 
a minimum delay angle for the converter, 

2,3,2 Firing Control Scheme with EPC 

An equidistant pulse firing scheme with pulse frequency 
control has been suggested to overcome the problem of 
harmonic instability which is inherent v^ith IPC [l6]. In the 
present context of simulation, only rectifier terminal is 
equipped with EPC, Fig. 2.8 shows a simplified block diagram 
of tlie controller from which the state equation can be 
written down, 

A constant slope ramp function is generated starting 
from zero at each firing instant. The ramp function is 
compared to the sum of the control signal (v^) and the bias 
voltage and pulse is initiated at the instant of equality. 

The bias voltage is a voltage proportional to T/6 in case of 
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6 pulse operation and T/l2 in case of 12 pulse operation, 

where T is the time period of ac voltage. When v =0, which 

c 

corresponds to the steady~state operation, the firing pulses 
are generated at every 60^ or 30*^ depending on 6 or 12 pulse 
operation respectively. To ensure the minimum delay angle 
operation for the rectifier the firing is inhibited until 
the commutation voltage has reached a defined minimum 
positive value. The details of this scheme are described in 
reference [l7], 

2.4 REPRESENTATION OF DC NETWORK 

The dc netv;ork consists of transmission line which is 
represented by a -n; model. The dc current of each converter, 
which is obtained in terms of state variables, is injected 
into the dc network at converter terminals. State equations 
for the dc network can be easily written dov;n and are time 
invariant. The voltage V in Fig, 2,1 is directly obtained 
as a capacitor voltage v,?hich is chosen as a state variable 
in the dc network equations. 

2 . 5 CONCLUSION 

In this chapter a model of a 12 pulse converter is 
developed for dynamic digital simulation. Both constant 
current control and constant extinction angle control are 
represented in detail. 
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CHAPTER 3 

SIMULATION OF HVDC LINK 

3.1 INTRODUCTION 

Based on the models of the converter, its associated 
controls and dc network, described in the previous chapter, a 
digital-computer program has been developed to simulate 
2 terminal HVDC link with 6 or 12 pulse converter at each 
terminal. Test simulations both under steady state and 
transient conditions have been carried out primarily to 
demonstrate the program capability. Commutation failure 
at the inverter terminal has also been simulated, 

3.2 COMPUTER PROGRAM 

Fig, 3,1 shows the flow-chart of the program which is 
written in FORTRAN and needs minimum data preparation. It 
has a modular structure which facilitates further program . 
augmentation. 

The dynamic equations of the entire system except those 
of the converter are time invariant. The time varying nature 
of these converter equations is due to the changes in the 
converter conduction pattern which occurs whenever a valve 
ceases or begins conduction. The two conditions which 
have to be satisfied, when a valve begins conduction are s 




Fig* 3»1 Flow chaa^t of the program 
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i) presence of gate pulse (which is the output of the 
controller routine), and 

ii) positive voltage across the valve, 

A valve ceases conduction if the valve current falls below 
zero and the exact instant of this extinction is computed 
by linear interpolation using valve current measurements. 
Each time the converter changes its status, equivalent 
circuit constants are evaluated (refer Chapter 2). The 
voltage source in the equivalent circuit is a function of 
the converter ac bus voltage, and hence it is computed for 
each time interval. Since the equivalent circuit constants 
are dependent on the converter conduction status, the latter 
is checked at the beginning. of each integration time step. 

The system differential equations are solved using 
Modified Euler method. 

The salient features of the program are as follows ; 

i) Both 6 pulse and 12 pulse operation of the converter 
can be simulated, 

ii) In its present form the program incorporate three 
methods of converter control which are i 

a) constant a 

b) constant current control 

c) constant extinction angle control. 
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The firing angle control is based on either IPC or EPC 
scheme (see Chapter 2). However simulation, that has been 
carried out , 5 the inverter operation is assumed only under 
IPC firing scheme (as this does not require feedback of 
extinction angle). 

iii) The voltage across the valve, following the cessation 
of its conduction is monitored fox a period equal to 
the turn-off time. If during this period the valve 
gets forward biased it is put back into the conducting 
state, 

3.3 TEST SIMULATION 

A 2 terminal HVDC system, shown in Fig. 3.2, has been 
considered for simulation, to demonstrate the program capa- 
bility, In this sample system, the rectifier (terminal l) 
and the inverter (terminal 2) normally operate on constant 
current and constant extinction angle controls rospectively. 
Various test simulations are conducted to investigate the 
system performance both under steady state and transient 
conditions following the disturbance. Details of the system 
parameters and operating conditions are given in the 
Appendix A. 

3,3,1 Steady State Response. 

The steady state waveforms for the 2 terminal system 
using 6 pulse converter at each terminal (operating 



3Z 
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condition A) and 12 pulse converter at each terminal 
(operating condition B) are shovm in the Figs, 3.3 and 3.4 
respectively. In either case IPC firing scheme is assumed 
at both the terminals. Since this study mainly illustrates 
the program capability, optimization of controller parameters 
has not been attempted, Vi/hile the dc current is almost smooth 
for the 12 pulse convorterj there is some 6th harmonic 
ripple for the 6 pulse converter, 

3,3,2 Transient Response 

Figs, 3,5 and 3,6 show the commutation failure in case 
of 12 pulse converter due to a 5 percent dip in ac voltage 
of phase B at the terminal. The dip is initiated at the start 
of cycle 1 and normal condition is restored at the start of 
cycle 2, IPC firing scheme is assumod here (operating 
condition B) , 


The transient system responses have been studied for a 
12 pulse 2 terminal system. Following disturbances were 
considered. Change in current reference setting (with 
operating condition B) using both (A) IPC and (B) EPC 
firing schemes. 


i) 

ii) 


at both the terminals, from p.u,? 

^ref^ "" ^ref^ == ^ref2 " 

at both the terminals from = 1.2 p.u. i 



34 



cyclefi 






Fig. 3.6 SHDV/S THE COMMUTATION FAILURE IN THE BRIDGES (12 PULSE CONVERTED 
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The system response for the above cases are shown in 
Figs* 3,7 to 3,10, All the above figures show the 
variation of average do current and dc voltages over 
50 cycles of simulation. The averaging is done for every 
30 degrees. In all the cases, thf disturbance is initiated 
at the beginning of 15th cycle of the simulation. 

Discussion 

The responses illustrate the controller action 
following the step change in the current reference setting. 
For case (i), while vdth IPC firing scheme the oscillations 
are observed in dc current and voltage, with EPC firing 
scheme these oscillations are negligible. With the increase 
in the magnitude of the step change in the current reference 
setting (case ii ) the oscillations in the dc current and 
voltage are noticeable with both IPC and EPC firing schemes. 
In the cases studied EPC appears to be more effective than 
IPC firing scheme. 

3.4 CONCLUSIOM 

In this chapter, a computer program developed for the 
dynamic digital simulation of HVDC system is described. The 
program employs both the individual phase control and egui— 
distant pulse control firing schemes for converter control. 
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Results of the various test slr»ulatf„ , 

m/nr 4 . s^ihlations for 2 terminal 

WDC system are presentoH r,,,- . 

prooented primarily to Illustrate the 

capahility of the program and they clearly show the need 

.*ose performance is 

influencGd by tho <;vc+cim 

y uie system parameters. 
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CHAPTER 4 

SIMULATION OF HVDC LINK WITH AC REPRESENTATION 

4.1 INTRODUCTION 

DGtailod representation of the ac system and harmonic 
filters is required to investigate the effect of interaction 
between ac and dc system dynamics on the system behaviour. 

The performance of the iiVDC system, particularly under 
transient conditions dotoriorates drastically when the ac 
system short circuit ratio (SCR) at the converter bus is low. 

A state space model of the ac system alongwith the 
harmonic filters has been formulated. Both steady state and 
transient studies have been carried out with detailed ac 
representation, after suitably augmenting the previous 
program described in the previous chapter, 

4.2 AC SYSTEM MODEL 

The schematic 39 representation of the ac system and 
harmonic filters feeding a particular converter is shown in 
Fig. 4.1. The ac system is represented by an ideal source 
voltage (Ej behind a T equivalent circuit for each phase 
of the ac network. The inductance (L^) and resistance (Rg) 
at the fundamental frequency are determined from the 
knowledge of SCR at the converter ac terminal and impedance 
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angle. Rg represents the effect of damping due to loads 
within the ac system. The shunt capacitance (C ) is 

3 

de tennined from the reactive power reguirement of the 
converter. The effect of the converter is represented as 
current sources feeding into the ac network 

The current through inductors in the links and voltages 
across capacitors in the tree of the graph are chosen 
as the state variables to define the system state at any 
given time. The generalised state equation of the ac 
system will be of the form 
• 

X = Ax + Bu (4,1) 

where u is the input, x is the state vector and A and B 
are constant mintrices. The ac system state and output 
equations are given in the Appendix B. 

4,3 INTERFACING THE CONVERTER AND AC SYSTEM 

Power from the ac bus is fed to the bridges of the 
12-pulse converter through two transformers each having an 
off-nominal tap sotting (a). One transformer has wye- 
wye connection vrtiero as the other has a delta-wye connection 
to obtain the desired 30 degrees phase difference between the 
two sots of voltages. The voltages and currents on the two 
sides of the wyo-wye connected transformer are related as 
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'2 ®bs ^ ®®3 ~ ® 


cs 


iY 


ap 


ai^Y j = ai Y 

9 » -*-j.p ai^ 


and the corresponding relations for the delta-wye connected 
transformer are. 


A 

'as 


f§ - ®2) 


= T§ (®2 - 


®cs - P (®3 “ ®l) 


■I ^ — - ( O’? ^ j. •4^\ 

^ap ■" Ji '■'^^9 ^7^ 


i6 


icp == yI ^^^7 + ^9^ 


The schematic diagram illustrating the interfacing of a 
converter and the ac system, is shown in Fig, 4.2. Here, 
and Ig are defined as follows 


1 1 * 4 « 1 ^ 

■^A ”* ap ap 

I =3 iY + i ^ 

•^B cp ^ cp 

The ac voltages ej^, eY^, e^^, e^J^, e^^ and e^^^ are 
used in the converter model (refer Fig. 2,1). The ac 
currents (i^, ig and ig) are obtained from the converter 




FIG. 4.2: INTERFACE BETVJEEN AC SYSTB4 AND CONVEFITER ]V10DEL 
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model. Since the program (described in the previous 
chapter) is developed on modular approach, it is anaugmented 
to include the ac system equations also. From the knowledge 
of the source currents (I^ and Ig) , equations B.l toB.7 
are solved at each integration time step to determine the 
voltages e^, e^ and 63 (i.e., e^ , refer Appendix B). The 
This in turn establishes transformer secondary voltages 
(Fig, 4,2) which are used subsequently for the solution of 
the converter dynamic equations. The- transformer tap is 
assumed fixed at tlie initial operating value determined 
from the load flow calculations, 

4.4 TEST SIMULATION 

The program capability is illustrated through simulation 
of a single terminal and two terminal HVDC system shown 
respectively in Figs, 4,3 and 3,2. In the AC system confi- 
guration shown in Fig, 4,1, each of the block marked 'F* 
contains 4 series tuned filters (i.e, 5th, 7th, 11 th, and 
I 3 th harmonic filters) and a second order high pass filter 
(refer Figs, B ,1 and B .2 for filter configuration). The 
detail system parameters and operating conditions are given 
in Appendix c. 

4,4,1 Steady State Response 

The steady state waveforms for a single converter, 
acting as a rectifier, with ac system of SCR equal to 15 and 
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i-ig* 4.3 SIMCLE 'fERIwIMAL MVDC SYSTEM 
FEEDING A CONSTANT VOLTAGE 
SaJilCE 





with filters included (operating conditions A) assuming 
IPC and EPC firing schemes, are shown in Figs. 4.4 and 4,5 
respectively. From Fig, 4.4 it is observed that the firing 
instants are not equally spaced, but not so for the case of 
EPC firing scheme shown in Fig, 4,5, For the same system 
without filters (operating conditions B) the steady state 
v/aveforms assuming IPC and EPC schemes are shown in Figs, 

4,6 and 4,7 respectively. The difference in steady state 
waveforms, with and without filters, is quite significant 
for IPC firing scheme but not for EPC scheme. Thus, it is 
evident, due to unequally spaced firing instants, IPC has an 
inherent tendency to enhance harmonic generation, 

4.4.2 Effect of Harmonic Filters 

The effect of harmonic filters can be seen from the 
test simulation results given in the Figs, 4,8 to 4,11, The 
difference in the ac currents and ac voltage of phase A with 
and without filters assuming IPC firing scheme are shown in 
Figs, 4.8 and 4,9 respectively. The above waveforms with 
EPC firing scheme, are shown in Figs. 4,10 and 4,11, Vtfhen EPC 
is used, hardly any disparity can be observed between the ac 
injocted currents- , Ip^ with and without filters. The 

above Figs. 4,10 and 4,11 also depict the nature of the 
injected current, CENi .u ,., ' 
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FIG. 4.4 STEADY STATE WAVEFOfiMS FOR SIKGLE 

TERMINAL SYSTEM WITH AC SYSTEM{ SCBr»15) 
AND FILTERS 

IFC FIRING SOHEJ'lE IS EMPLOYED 
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Discussion 

The steady state responses, for a single converter 
system, show that the filters r.re necessary if I?C 
firing scheme is employed at the rectifier. Thus, the 
distortion in the ac voltage has least effect on EPC firing 
scheme, whereas IPC is more prone to introduce harmonics as 
its firing pulse generation is based on the zero crossing 
of the commutation voltages. Comparison of the responses 
given in Figs. 4.12 and 4.13 (assuming IPC and EPC schemes 
respectively) shows that IPC firing scheme introduces low 
frequency oscillation in the dc voltages but EPC scheme does 
not, 

4,4,3 Transient Response 

Transient responses are determined for a 2 terminal 
HVDC system where ac system representation is considered 
only at the inverter terminal. The disturbance chosen is 
the change in current reference setting to observe the 
system responses v/ith IPjC and EPC firing schemes when the 
inverter is feeding a strong and a weak ac system. 

For 2 terminal HVDC system 

Aj with IPC firing schemes at both the terminals 

i) SCR = 15 (operating condition C) 

ii) SCR = 7 , (operating condition D) 
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B: with EPC at rectifier and IPC at inverter 

i) SCR = 15 (operating condition C) 

ii) SCR = 7 (operating condition D) 

In both the cases the change in current reference setting 
^ref from 1,0 to 0,75 for terminal 1 and from 0,8 to 
0*55 for terminal 2, The system responses for the above case; 
are shown in Figs, 4.14 to 4.17. For all the cases the 
disturbance is initiated at the 10th cycle of the simulation. 

Discussion 

The transient responses for a two terminal system show 
that, follov;ing a change in the current reference setting, 
oscillations in dc current and voltage are present v;ith IPC 
firing scheme at all terminals. It becomes quite significant 
if the converter is connected to ac system of low SCR (7), 

This undesirable influence of the ac system and filter 
parameters may perhaps be eliminated by proper controller 
design and use of control filters, V^ith EPC at the rectifier 
alone, oscillation are found to reduce considerably and 
current settles down to the new sotting in a reasonably 
short time, 

4.5 CONCLUSION 

In this chapter a state space model of the ac system 
including harmonic filters is described. This is interfaced 
with the converter model described in Chapter 2, The computer 
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program developed earlier is suitably changed to incorporat 
the ac system representation. Steady state and transient 
responses of single and tv;o terminal system are presented 
respectively. The responses not only show the influence 
of ac systems and filters but also points out that dynamic 
performance of HVDC system connected to ac system of lov^ SCR 
can be improved greatly by employing EPC at the rectifier 


alone 
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CHAPTER 5 

CONCLUSION 

This thesis is mainly devoted to the development of a 
12 pulse converter model for dynamic digital simulation of 
HVDC system to investigate the control system response. All 
the subsystems including the dc network, ac system with 
filters, converter and the control system have been 
represented in detail. 

The converter is represented by an equivalent circuit of 
a variable voltage source behind a variable impedance. Both 
these quantities are functions of the valve conduction pattern^ 
This equivalent circuit is derived based on graph theoretic 
analysis. The model facilitates simulation of all possible 
modes of converter operations having continuous dc link 
current. The converter model is general enough to handle non- 
ideal valve characteristics and unsymmetrical source impedances; 

' i 

The computer program developed incorporating the proposed 
converter model is modular in structure. As a result, it 
increases flexibility in the representation of individual 
subsystems. 
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The digital simulation was carried out on DEC 1090 
system and to simulate one cycle of ac waveform the program 
takes 4 to 12 seconds of CPU time depending on the system 
details. This time can be further reduced by program 
optimization. Due to lack of standard test data, the results 
obtained in the thesis could not be compared witii the results 
obtained from other methods of simulation. 

Conclusions drawn from the test simulation results 
obtained in Chapters 3 and 4 are as follows s 

1, EPC is more effective than IPC particularly when 
operating in a system without ac filters or when the 
converters are connected to a weak ac system, 

2, The system performance is critically dependent on the 
controller configuration and parameters. 

The major thrust in this thesis was on the development 
of the necessary software for the dynamic digital simulation 
of 12 pulse converters. However, the following aspects have 
not been taken into account : 

1, Detailed representation of synchronous generator and 
transformer with their associated magnetic non- 


linearities 
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2, The study of the performance of static var compensation 
in the integrated ac/dc power system, 

3, Program optimization using more efficient numerical 
technique. 
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APPENDIX A 


The dc base quantities are chosen as 

DC base voltage is ideal no-load dc voltage of 
a bridge = 100 KV 

DC base current = 1 KA 

DC base impedance = 100 Ohms 

The system parameters are adapted from reference [ j ] 

Two terminal system 

Rated Power of the DC link = 240 MIV 
AC system frequency = 50 cycles 


Bridge transformer s 


Resistance R^^ = R ^2 = 0,5 ohms 
Commutating Reactance 


= 6,283 Ohms 


Smoothing Reactor s 


Resistance 

Reactance 


^dl " ^d2 
^dl ^d2 


=0.1 Ohms 
= 314.159 Ohms 


Transmission line : 


Total resistance = 8,64 Ohms 
Total Inductance = 0,50148 H 
Total Capacitance - 54,1645 iiF 
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Controllers ; 

a) Individual phase control 

~ K 2 = 2070 Ohms units 
Tci = =0.1 seconds 

b) Equidistant phase control 

= 0.001 rad/amp-sec. 

T 1 /T 2 = 2.0, = 0.0033 seconds 

RTvlS AC voltages (L-N) ; 


= 42.78 KV, = 38.51 KV 

Valve turn off time = 3^ 

Operating Conditions ; 

A. = 1^12 ~ 1000 Amps, 

%in=^°» Tc = 10" 

= 26.22*^, = 148.39° 

B. ~ ^d2 ~ 1200 Amps 

a . = 5° , Y = 10° 

min * 'c 

= 26.11° y = 145.59° 


Terminal 1 (rectifier) on constant current control and 
Terminal 2 (inverter) on CEA controls. 
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appendix b 


From Fig. 4.1 it can be found out that three atate 

variables are necessary and j. ^ 

^ to describe the 

ate of each phase (except filt 0 r<i^ ana 4 -k ai. 

^ ^ J-tcTs) and the three genera- 

lisea state equations for jth phase are 


Ij 


■2j 


where 


P 






if 




■IJ - s 
J = 1 

j = 2 


(B.l) 

(B.2) 
(B. 3) 


I = T 
P, B 


IX 


j = 3 


~ ^^stantanc. 


n 


■Fj 


where i^ 


2 : 

K=1 


'■ous ideal source voltage of phase j, 
j = 1,2, and 3, 


K / 

j +(ihf j + (ej - '^hfj ) /Ry^f ) 


tfj IS the current through the Kth tuned flits 
is current through the high-pass filter. 


and 
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The block marked *F' in Fig, 4,1 represents a set of 
series tuned filters (5th, 7th, llth and 13th) and a high 
pass second-order filter, A series tuned filter as shown 
in Fig, B,1 can be represented by two state equations as 


^tfj ~ “ ^tf Hfj^/^tf ■ 

^tfj ^tfj/^tf 

while the state equations describing the high-pass filter 
shown in Fig, B,2 can be Vi/ritten as 

ihfj = [Sj - ''hfj^/Lhf (B.6) 

^hfj = C^hfj tB.7) 


K K 

and are the capacitance, inductance and 
resistance of the Kth harmonic tuned filter, and and 

are the capacitance, inductance and resistance of the 
high pass filter, 

els» the converter ac bus voltages are the output of 
the system whereas I^, Ig and are the input. 
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FIG. B.l: SERIES lUNED FILTER 




e . 
3 


9 



Seco 

FIG. B.2: SECOND-ORDER HIGH PASS FILTER 
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APPENDIX C 

All the ac system parameters are given with reference 
to the Fig, 4.1, 

Single Terminal (Rectifier) 

Transformer : 

Tap a = 1.0 

Resistance R^ = 0.5 Ohms 
c 

Reactance X = 6,283 Ohms 
Converter : 

Valve turn off time = 3^ 

AC Bus Voltage e. = 42,78 KV 

rT-.w> J 



DC current = lOOO Amps 

Constant Voltage Source V = 166,449 KV 

c 

®min ~ * ^c =10.0 

= 26.9^ 

Smoothing reactor 
Resistance R^ = 0,1 Ohms 
Reactance X^ = 314,159 Ohms 
Controller i 

a) Individual phase control 
Kl = 20^0 Ohms 
T^ = 0.002 seconds 

Q 
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b) Equidistant Pulse Control 
Kj_- C.CX}1 Xiid/arap-sec 
T^/T 2 = 2 f T 2 = 0,0033 seconds. 
Operating Condition ^ 

A ; With filters ; 

AC System j 

Short circuit ratio SCR = 15 
Impedance angle = 85° 

Resistance R^ = 4,16585 Ohms 

Reactance : 

Inductive = 0,74044 Ohms 

Capacitive = 263.676 Ohms 

Ideal ac source voltage 
Ej = 43.88 KV 

©j = 13.20948° 

B : Without filters 
AC Sys tern : 

Short circuit ratio SCR = 15 

Impedance angle = 85° 

Resistance R = 4,3072 Ohms 

Reactance i 

Inductive = 0,76559 Ohms 

s 

Capacitive C = 263,676 Ohms 
Ideal ac source voltage 
Ej s 44,67 KV 

Qj « 13,17981° 
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Two terminal j 

AC system is represented at Inverter only. 

Only AC system and operating conditions are given below. 
Other system parameters are same as given in Appendix B. 
RMS ac source voltage (L-N) at rectifier 
e^ = 42.78 KY 

0 ^ = 0 . 0 ° 

at inverter 

e^ = 38.51 KV 

02 = 0 . 0 ° 

DC current == I.O KV 

== 26.9° , = 148.39° 

Operating conditions : 

C : Short circuit ratio SCR = 15 

Impedance angle s= 85° 

Resistance R = 3.04 Ohms 

s 

Reactance ; 


Inductive = 0.605 Ohms 
Capacitive C^ = 263,676 Ohms 
Ideal AC source at invertor 


Ej = 39.088 KV 

04 = 7.03597° 


Transformer tap 




D ; Short circuit ratio SCR = 7,0 
Impedance angle =s 85^ 

Resistance R^ =: 6.1859 Ohms 
Reactance t 

Inductive = 1.0995 Ohms 
Capacitive Cg= 263.676 Ohms 
Ideal AC source at inverter 
Ej = 39.112 KV 
= 6.004^ 

Transformer tap a, = = 1.0 

Filter Parameters 


Those parameters arc taken from the roforonco [.,] and used 


for all the 

test simulations 

«k 


Harmonic 

No. 

Resistance 

Ohms 

Inductance 

mH 

Capacitance 

(iF 

5 

1.11413 

362.2 

11.1346 

7 

0.801103 

181.6 

11,13748 

11 

0.504805 

726,5 

11.5.102 

13 

0,4225 

520.4 

1 1.5331 

H.P. 

16.1428 

3,0041 

1 1 . 57' 5 1 



